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The dynamics of hexagon patterns in rotating systems are investigated within the framework of 
modified Swift-Hohenberg equations that can be considered as simple models for rotating convection 
with broken up-down symmetry, e.g. non-Boussinesq Rayleigh-Benard or Marangoni convection. In 
the weakly nonlinear regime a linear stability analysis of the hexagons reveals long- and short-wave 
instabilities, which can be steady or oscillatory. The oscillatory short-wave instabilities can lead 
to stable hexagon patterns that are periodically modulated in space and time, or to a state of 
spatio-temporal chaos with a Fourier spectrum that precesses on average in time. The chaotic state 
can exhibit bistability with the steady hexagon pattern. There exist regimes in which the steady 
ON ■ hexagon patterns are unstable at all wavenumbers. 
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I. INTRODUCTION 



Patterns arising from hydro-dynamic instabilities have become a paradigm for the investigation of dynamical sys- 
tems with many degrees of freedom. In particular Rayleigh-Benard convection has been studied experimentally and 
theoretically in great detail. In recent years the focus of investigation has turned to patterns that are spatially and 
possibly also temporally complex. Various types of such states of spatio-temporal chaos have been identified. Great 
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I interest has, for example, found the spiral-defect chaos that is observed in convection of fluids with not too large 
. Prandtl number |^,|3|]. It is characterized by the appearance and disappearance of spirals in the convection patterns 
that are separated by domains of straight or slightly curved convection rolls. Spiral-defect chaos has not been found 
G\ . directly at onset of convection. Thus, it is not accessible by weakly nonlinear theory. From a theoretical point of view 
another type of spatio-temporal chaos has therefore been very attractive since it arises directly at onset. It is observed 
in rotating convection and is due to the Kiippers-Lortz instability through which all straight-roll states become unsta- 
ble to rolls of a different orientation if the rotation rate is sufficiently large [Q . In a simple three- mode model for rolls 
rotated by 120° with respect to each other this instability leads to a structurally stable heteroclinic orbit connecting 
the three sets of rolls |^||. In large aspect-ratio experiments @-||] and numerical simulations of Ginzburg-Landau 
equations JhJ and of order-parameter models |ll|]l2f | the switching between the different orientations looses coherence 
across the system and domains of almost straight rolls of different orientations form that persistently invade each 
other in an irregular manner. 

The character of a spatio-temporally chaotic state is expected to depend strongly on the planform of the underlying 
pattern. Thus, spatio-temporal chaos arising from instabilities of stripes (or rolls) will differ from that of square or 
hexagonal patterns. For instance, in hexagon (or square) patterns the bending of rolls that is widely observed in 
Rayleigh-Benard convection is strongly suppressed. Instead, the appearance of disorder in a hexagonal pattern might 
resemble more the melting of a 'crystal'. There have been only a few studies of transitions to spatio-temporal chaos 
from square or hexagon patterns. Experimental investigations have been made on parametrically driven surface waves 
p3| , pd| and on a chemical system [|l5] . 

In the present work we consider hexagon patterns in systems with broken up-down symmetry and study the effect of 
a breaking of the chiral symmetry on their stability and dynamics. It is motivated by the complex dynamics ensuing 
from the Kiippers-Lortz instability of convection rolls, where the chiral symmetry is broken by a global rotation of 
the system. From a general weakly nonlinear analysis of the three modes making up the hexagon patterns it is known 
that a breaking of the up-down symmetry transforms the structurally stable heteroclinic cycle that results from the 
Kiippers-Lortz instability within the three-mode model into a periodic orbit of oscillating hexagons f]Jl§|l3 . This 
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periodic orbit arises from a Hopf bifurcation, which replaces the steady instability of the hexagonal pattern that leads 
to rolls D@. 

We are interested in the question which additional instabilities can arise if one goes beyond the three-mode model 
and allows arbitrary perturbation modes. Of particular interest are situations in which such perturbations lead 
to persistent dynamics. One possible approach is to extend the three-mode amplitude e qua tions to three coupled 
Ginzburg-Landau equations that also allow for a slow spatial variation of the amplitudes . The slowness of the 
space dependence implies, however, that all perturbations to the hexagons are restricted to wavevectors close to 
those of the hexagon modes. In particular, only modes making a small angle with respect to the hexagon modes can 
be treated. Furthermore, the Ginzburg-Landau equations cannot be expected to be suitable for general numerical 
simulations of unstable hexagon patterns in large systems since the choice of the three carrier modes making up the 
hexagons strongly breaks the isotropy of the system. Only in the unlikely case that the perturbation wavevectors 
remain close to the wavevectors of the hexagon modes for all times would the Ginzburg-Landau equations be adequate. 

The goal of our investigation is to explore what type of instabilities can appear due to the breaking of the chiral 
symmetry and whether they can lead to persistent regular or irregular dynamics. It is therefore important to allow 
arbitrary perturbation modes and to preserve the isotropy of the system. Thus, we investigate a simple model of the 
Swift-Hohenberg type jlq^Oll . Such models have been widely used to address general questions regarding the dynamics 
of patterns (e.g. [|1T Jl9|jl 2j , |2lf p3|). If the patterns arise from a long- wave instability equations of the Swift-Hohenberg 
type can be derived systematically in a long-wave analysis. This has been done, for instance, for buoyancy-driven 
convection with asymmetric boundary conditions and for Marangoni convection J25|,[32| . In both cases the onset 
of convection occurs at long wavelengths if the boundaries are poor conductors. If the wavenumber at onset is finite 
the Swift-Hohenberg cannot be derived rigorously from the basic equations (e.g. Navier-Stokes equations), but they 
can be obtained as approximate order-parameter models by expanding the nonlinear interaction terms (e.g. p9|j20[ |). 

We perform a general linear stability analysis of the weakly nonlinear hexagon pattern and pay particular attention 
to destabilizing modes that are rotated with respect to the hexagon pattern. In the absence of rotation such general 
stability analyses of hexagonal patterns have been performed for the full fluid equations for Marangoni convection 
p6| as well as within coupled Ginzburg-Landau equations p7| . Some results are also available for weakly nonlinear 
Marangoni convection in the presence of rotation p8[ . We complement the stability analysis with direct simulations 
of the order-parameter equation to investigate the nonlinear behavior resulting from the instabilities. We find that the 
hexagonal pattern can become unstable at all wavelengths, resulting in spatio-temporally chaotic states. Moreover, 
even for parameters for which stable hexagon patterns exist, persistent periodic and irregular dynamics are found, 
implying bistability of ordered and disordered states. 

The paper is organized as follows. In Section II we introduce two order-parameter models and set up the linear 
stability analysis of the weakly nonlinear hexagons. In Section III the results of that analysis are presented. We 
discuss the different types of instabilities and the resulting stability regions of hexagonal patterns. In Section IV 
numerical simulations of the order-parameter equation are used to confirm the stability analysis and to investigate 
the nonlinear evolution ensuing from it. Conclusions are found in Section V. 



II. THE MODEL AND LINEAR STABILITY ANALYSIS 



Spatially periodic hexagon patterns with small amplitude arising from a weakly transcritical bifurcation can be 
described systematically by coupled equations for the complex amplitudes A, B, and C of the three stripe (roll) 
components 

|^ - R 2 A + 2CB*C* + (MAI 2 + f 2 \B\ 2 + f 3 \C\ 2 ) A, (1) 

where t 2 is a slow time, R 2 a reduced control parameter and the difference between the cubic coupling coefficients f 2 
and /a is proportional to the breaking of the chiral symmetry. The equations for B and C follow by cyclic permutation. 
To investigate the stability of these hexagons with respect to spatially varying perturbations the amplitudes could 
be allowed to be space-dependent which would lead to the inclusion of spatial derivatives in ([[]). The amplitude 
description breaks, however, the isotropy of the system and requires that the space dependence of the amplitudes be 
slow. This implies that only perturbations with wavevectors close to those making up the hexagon can be considered. 
In this study we are interested in particular in perturbations that are rotated with respect to the hexagons by an 
arbitrary amount and in dynamical, disordered states that may exhibit an essentially isotropic wavevector spectrum. 
Therefore the isotropy of the system needs to be preserved and Ginzburg-Landau-type amplitude equations are not 
sufficient to investigate the questions we are interested in. 
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For a first exploration of the impact of a chiral symmetry breaking on side-band perturbations of hexagons and 
their nonlinear evolution it is appropriate to preserve the isotropy of the system by considering a simple model 
equation rather than the full hydrodynamical problem. Detailed investigations of specific physical systems can then 
be performed as a second step guided by the results of the model equations. 

In a number of studies Swift-Hohenberg-type models have provided valuable insight into the dynamics of patterns 
(e.g. [|l7Ul9|Jl^ , pl| -^3|). We consider here a minimal model of that type for an order parameter ip(x,y,t). In order 
to obtain stable hexagons in the absence of rotation we introduce a quadratic term atp 2 that breaks the up-down 
symmetry ip —> —ip. The cubic scalar term ip of the original Swift-Hohenberg equation provides nonlinear saturation. 
To break the chiral symmetry nonlinear gradient terms have to be introduced. At leading order in the amplitude and 
the gradients this introduces an additional quadratic term and we are led to consider the model 

d t ijj = Rip~ (V 2 + l) 2 ip - ip 3 + atp 2 + je z -{Vipx V(y 2 iP)), (2) 

where e z is the unit vector perpendicular to the (x, ?/)-plane. In these equations the bifurcation parameter is R and the 
wavenumber of the bifurcating mode has been scaled to 1. The strength of the chiral symmetry-breaking is measured 
by 7. 

It should be noted that the quadratic gradient term in (|2[) has the somewhat non-generic feature that it generates 
harmonics only from the interaction of modes with wavevectors of different magnitude. As a consequence, in the 
small-amplitude regime (Q), for which a has to be assumed to be small, it does not modify the cubic coefficients fi 
in ([j]) and has therefore to this order no effect on strictly periodic hexagon patterns. In particular, the instability 
of hexagons that occurs with increasing R is not oscillatory and does not lead to the periodic orbit of oscillating 
hexagons identified in [fj]]l6|] (cf. Fig.|J). Instead, the instability remains steady and leads directly to the roll solution 
as is the case in the absence of rotation. In this paper we focus on the stability of the steady hexagons with respect 
to side-band perturbations below that instability. Thus, the existence of the secondary branch of oscillating hexagons 
is not of central importance. We present most results for the minimal model (||). 

The non-genericity of the quadratic gradient term in (|^) can be addressed by introducing additional terms. One 
possibility is to introduce a quadratic term, e.g. (3(\7ip) 2 : that allows a to be of 0(1) as long as a and (3 together lead to 
£ <C 1 in (Q). Alternatively, a cubic term that breaks the chiral symmetry can be introduced, e.g. 172 e 2 '(Vx (V 'tp) 2 \/ 'ip) . 
We thus have the extended model 

d t i> = Rip - (V 2 + l) 2 ^ - V 3 + aV 2 + /?(VV0 2 + 7 e z ■ (W> x V(V 2 tp)) + g 2 e z • (V x (V?/>) 2 V</>), (3) 

whose results we discuss briefly. 

While we are using (|^) and (|§|) as simple model equations, equations of that type can be obtained from the basic 



physical equation (e.g. the fluid equations) by expanding the interaction term in Fourier space |9|j30||31j]. In general, 
one then obtains additional nonlinear gradient terms. Similar equations can be derived systematically as long-wave 
equations (e.g. for buoyancy-driven convection with asymmetric boundary conditions j24j] or for surface-tension driven 
convection without [^5| and with rotation J32|). 

We perform a weakly nonlinear analysis of (^|) and (|^), and start with a small amplitude expansion of ip on a 
hexagonal lattice, 

if, = e ^Ae ikx + 5 e lfe (-z/2+yV3/2) + Ce ik(-x/2- y V3/2) j ^ 

1 ,2 f ?-) 1 n J.2kx 1 7-1 ik(—x+V3y) , rj „ik(— x — V3y) 

+ £ j ^000 + ^20oe + i^o2oe + ^oo2e 

1 n . „.t(3i;/2- s 4/2) , n ik(3x/2+ v V3/2) , jy ikV3y 1 

+ c.c. + h.o.t., 

where e<l and k is the wavenumber of the pattern. The Fourier amplitudes are functions of the slow time r 2 — t 2 t. 
The subscripts on the harmonics D\ rnn indicate how the Fourier mode is expressed in terms of the O(e) modes; the 
first, second, and third subscripts indicate how many times e lkx , e lk (- x / 2 +v^/ 2 ) j and e ik (- x / 2 -yV3/2) are f ac ^ ors of the 



Fourier mode, respectively. Bars over the subscripts indicate complex conjugates. For example, D 1Q -e 



ik(3x/2+yV3/2) 

^ .. ^ „ „ „ iJr ..„, ~ 1UA ~ 

contains one e lkx -mode and one complex conjugate of the e lk ^~ x l 2 ~ v ^ 3 ~l 2 ^ -mode. 

Inserting the expansion (Q) into the model equations (||) and (^|) and eliminating the terms involving Di mn at 
0(e 2 ),in each case one obtains at 0(e 3 ) a set of amplitude equations of the form (|l|) for the amplitudes A, B, and C. 
In the case of the minimal model (||), the coefficients of ([j]) have the values 

/i=3, / 2 = 6, / 3 = 6, (5) 
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and R2 = (R — (k 2 — l) 2 )/e 2 . Note that R is expanded about the neutral curve for the specified hexagon wavenumber 
k, and not about R = 0. To obtain a homogeneous scaling, we have assumed the bifurcation to be weakly transcritical, 
and hence 

C = a/e = 0(l). (6) 

Note that the quadratic rotation term 7^ • (Vi/> x V(V 2 f/0) does not enter the amplitude equations ([!]) if a <C 1. In 
terms of the original variables, the amplitude for the steady hexagon solution is to this order given by 



, r> r -a - V* 2 ~(R- (fc 2 - I) 2 ) (A + h + h) ~ 

(Jl + J2 + J3J 

In the case of the extended model (||), the coefficients of (0) have the values 

C= / Q+ 2*!\/ e = (l), (8) 



/l = 3 -{ 2 (6 fc2 !i5^) + F^F} /32fe4 ' (9) 



The effect of rotation is contained solely in f% and /3 such that 72(7,32) = 7, — 32)- Note that in (|9|,|10|) a has 
been replaced by /3 using (g). The hexagon solution is to this order given by 



+ ^) 2 - ( fl - ( fc2 - !) 2 ) (/1 + /2 + h) , s 

A = B = C = ih + h + h) ' 

In either case, to study the stability of the hexagon solution H(x,y) a general Floquet ansatz with modulation 
wavevector (q,p) is made, 

iP = H(x,y)+5e iiix+i ™P(x,y,t), (12) 

where P(x, y) has the periodicity of the hexagonal lattice and 6 <C 1. Being interested in the weakly nonlinear regime, 
we expand H(x, y) and P(x, y) in e and obtain 

% j ) = e ^Ae ikx + Be lk ^ x/2+y ^ /2 ' ) + Ce ik! -~ x/2 ~ v ^ /2 ' ) } (13) 
+ e 2 { A,oo + ^20oe l2fea; + D 02Q e ik (- x+VS ^ + D 00 2e tk ^ x ~^ y) 

+ D _ Jk(3x/2-yV3/2) , n ife(3a:/2+8/v^/2) , r) -ifcV^iA 

-r-L/ 110 e -r ^i i e > ^on e / 

+ e 2 5 |i7 - e i ((«+ fc / 2 ) :E +(P- fc \ / 3/ 2 )?') + i/ 100 e i((9+fc):E+TO) + H Q0 -e l{( - q+k/2)x+( - p+k ^ /2)v) 

+ H Tme i((l-k)x+py) + p owe i((q-k/2)x+(p+kV3/2)y) + p^ e i((q-k/2)x+(p-kV3/2)y) j 
+ C.C. + /l.O.t. 

with the amplitudes of the perturbation modes depending on time. Here the notation for the O(e) and 0(e 2 ) modes 
is as before. The subscripts on the 0(e 2 S) amplitudes indicate how many times the e^ qx+pv ^-mode was multiplied 
by the indicated Fourier components that make up the original hexagon solution to create the harmonic mode. The 
reason for labelling some of the 0(e 2 5) modes as H and some as P will become apparent below. In the expansion 
( |l3| ) we have combined the modulation wavevector (q,p) with the basic wavevector (k, 0) of the hexagon pattern into 
the perturbation wavevector (q = k + q,p = p). Thus, long-wave perturbations are characterized by (q,p) being close 
to (fc,0). 

As long as the wavevector (q,p) of the perturbation mode Pooo is not close to the wavevector of one of the hexagon 
modes the harmonics Hi mn and Pi mn of Pooo are bounded away from the critical circle (cf. Fig.[l] below). Thus, their 
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damping is 0(1) and they can be eliminated in favor of Pooo and the hexagon modes A, B, and C. Consequently, 
they are small (of 0(8e 2 )) compared to Pooo- This results in an evolution equation for Pooo alone. 

It turns out that the expansion ( |l3| ) yields unsatisfactory results for the growth rate of the perturbation as a 
function of 8. This is due to singularities that arise for 8 — and 8 = ±60°. For k = \J q 2 + p 2 = 1, as 8 — > 
the perturbation modes Poio and Pooi approach the critical circle and their linear damping goes to zero. Thus, 
they cannot be eliminated adiabatically any more as signified by singularities in the growth rate for Pooo- Similarly, 
when 9 — ► 60° the mode Polo approaches the critical circle. The influence of the singularities at 0° and 60° on the 
eigenvalues is so large that they become unreliable over the whole range 0° < 9 < 60°. 
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FIG. 1. Fourier space diagram showing the critical unit circle and location of the modes used in the final expansion of ip (cf. 
([[i])). Indicated with arrows is an example of the combination of an 0(e8) perturbation mode and an 0(e) original hexagon 
mode driving a 0(e 2 5) higher harmonic mode. 

FIG. 2. Fourier space diagram where 8 — > 60°. Here the modes H Q - , H j lt H002, and Pyoo approach the critical unit circle 
despite being of 0(e 2 S), and thus drive a singularity when solved for. 



In order to avoid the singularities, Poio and P001 are allowed to be of 0(eS), i.e. of the same order as the original 
perturbation mode Pooo- Then at 0(e 2 S) they drive the modes Pj 01 and Pj 10 leading to singularities in these modes 
for 9 — > 0. Hence, modes Pj Q1 and Pj 10 must also be allowed to be of 0(e5). In fact, all modes that come close to 
the critical circle as e l ^ qx+pv ^ approaches e lx must be kept at O(eS). They are labelled P; mn . This finally leads to the 
following expansion for tp (cf. Fig. Q), 

^ _ e ^Ae lkx + Be lk{ ~ x/2+v ^ /2) + c e * fe (-*/2-y\/3/2) j ^ 

+ eS | P aoe i{qx+py) + P 010 e i{ - {q ~ k/2)x+i - p+k ^ /2)y) + p Q0l e l ^ k / 2 y x+ ^ k ^/ 2 ^ 

i p_ p i((q~3k/2)x+(p+kV3/2)y) , p_ i((q-3k/2)x+(p-kV3/2)y) , p_ i((q-2k)x+py)\ 
T +110^ T -Tioi^ ' ^ 200^ J 

1 ,2 J n in „i2kx 1 n ik(-x+^/3y) . n ik(—x-V3y) 
+ e |^000 + ^200 e + ^020e + Dqq2& 

4. n _ p tk(3x/2-yV3/2) , n ik(3x/2+ V V3/2) , n _ p ikV3y\ 
' ^110 e ' ^101 e ' ^011^ J 

+ e ^5 ^H Jao e t{ ^ k)x+pv) + H 10Q e^ q+k)x+pv) + H oaJ e^ q+k/2)x+(p+k ^ 2 ^ 

-|- JJ _.^ e i(lx+(p-kV3)v) _|_ U QQ2 pi{{q~k)x+(p-kV3)y) _|_ ^_^ e i(( 9 -2fc):r+(p-fc\/3)?/) 

1 H- p i((q~5k/2)x+(p-kV3/2)y) , rr_ i((g-3fe)a+ TO ) , u_ i((«-5fc/2)z+(}H-*V3/2)y) 
' -"201° ^-"300^ '-"210^ 

+ iJ- 20 e i ^ 9_2fc ' ;r+ ^ +fev ^^' + # 020 g*((?- fc ) ;c +(j , + fe v / 3)3/) _|_ jj & i(.Qx+(p+kVS)v) 
+ ff^^^)} + c.c. + fc.o.i. 

Certain values of (<?,p) drive a resonance that is still not accounted for in the scaling of the modes. For k = 1, 
as (9 approaches 60° the expansion breaks down again because the modes corresponding to H n - n , P n - 1 , H- nri > and 
P002 are now near the critical circle (cf. Fig. ||), but are still kept at 0(e 2 S) in the expansion (|l4|). Similarly, the 
expansion breaks down as 8 approaches —60°. Due to the 60°-periodicity of the expansion, we may restrict ourselves 
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to 9 between —30° and 30°, and thus avoid having to do an altered calculation for these situations. If k = O.i 
a resonance occurs already for 9 = ±30°. Since k — .888 is smaller than any hexagon wavenumber we consider, 
restricting 9 to the range —30° to 30° is adequate. 

To summarize, in the expansion (|T^ ) all modes are retained that are excited at 0(e 2 ) or 0(e 2 5) by the perturbation 
modes or the hexagon modes. This is sufficient to determine the growth rates of the perturbation modes at 0(e 3 8). 
The expansion remains well ordered as long as the hexagon wavenumber and the perturbation wavenumber are close 
to 1 and the angle 9 is between —30° and +30°. 

We now substitute ([l4]) into (0), expand, and use equation ffih to find the perturbation growth rates for the minimal 



model. For the extended model (|3|) and (11) are used. To obtain a unified description for the whole range of angles 
we reconstitute the perturbation equations at 0(e 2 5) and 0(e 3 S) by substituting in for all scaled amplitudes, e.g. 
A = A/e, Pimn — Vimn/^S, etc., as well as the scaled variables, C = ot/e or £ = (a + /3/2)/e (depending on the model) 
and i?2 = (R — (k 2 — l) 2 )/e 2 , and the scaled times, dx — dt + ed Tl + e 2 d T2 . This leads to 

d T Vovo = (R~ (q 2 +P 2 - l) 2 )^ooo (15) 
+ B'Poio/oiofep, k) + C*V Qifoai{q,p, k) + ^ 9imn{q,P, k)\A\ 2 Vi mn , 

hnn 

where the sum goes over all six perturbation modes Vimn- The equations for the other amplitudes Vimn follow by 
cyclic permutation within the circle of modes (-Pood -Poicb Piw? -F20O) ^01 > -fboi)- It turns out that despite the 
assumption of weak transcriticality, @ or (||), the coefficients of the quadratic terms involving B* and C* are not 
small in general. While the coefficients go to zero for 9 — > and ( — 0, they are of 0(1) for finite angles 9. For finite 
angles and k close to 1, however, the modes Poio ancl ^001 to which T-ooo is coupled are small rendering the quadratic 
terms again of the same order as the cubic terms. 



III. LINEAR STABILITY RESULTS 



The general stability analysis (|15|) leads to 6 eigenvalues. To provide a context we first discuss the special case that 
the perturbation modes coincide with the original hexagon modes (i.e. (q,p) = (k, 0)). The analysis is then equivalent 
to that of the amplitude equations (|l|) , which can be rewritten as equations for three real amplitudes and three phases. 
Due to translation symmetry in two directions, two of the phases decouple and the eigenvalues corresponding to the 
translation modes are identically zero. All of the relevant behavior is contained in the three real amplitude equations 
and one phase equation. As long as the up-down symmetry is broken sufficiently strongly the remaining phase relaxes 
quickly to a stable fixed point, and only three real amplitude equations are needed, corresponding to three eigenvalues 
instead of six. This system has been thoroughly studied by Swift and Soward jl],[l6| . In Fig.|| the resulting bifurcation 
diagrams are sketched for the case without and with rotation. The new feature introduced by the rotation is a limit 
cycle of oscillating hexagons that branches off the hexagons and connects it with the general solution in which all 
three amplitudes have different magnitudes. The limit cycle arises in a Hopf bifurcation at 

p 8e 2 C 2 (/2 + / 3 + 4/i) tu2 a 

R Ho Pf = . „ +{k - 1) , (16) 

1/2 + J3 - 6Jl) 

with C given by (||) and (@), respectively. At Rnopf the real parts of two of the three amplitude eigenvalues become 
positive. Without rotation the eigenvalues are real and the instability leads through a transcritical, steady bifurcation 
directly to the unstable general solution and renders the rolls the only stable state. Note that because f 2 and fa 
only appear in the combination /a + .^3, Rnopf depends on neither 7 nor 172- The same is true of the saddle- node 
bifurcation at which the hexagons come into existence. In fact, 7 and g 2 are relevant solely for the imaginary parts 
of the eigenvalues. 



G 




FIG. 3. Bifurcation Diagrams (cf. Q). The hexagon solution gains stability through a saddle node bifurcation. Without 
rotation the hexagons become unstable through a transcritical bifurcation involving a general solution, where the three ampli- 
tudes are all nonzero but not all equal. With rotation the hexagons become unstable to a limit cycle through a Hopf bifurcation. 
The limit cycle in turn disappears through a collision with the general solution. 

In this paper we are interested in perturbation modes that differ in their wavevector from the hexagon modes (i.e. 
(q,p) ^ (k, 0)) and consider the dependence of the eigenvalues on (q,p). Thus, we deal with six branches of eigenvalues. 
Away from the saddle-node bifurcation and for sufficiently strongly broken up-down symmetry only four branches are 
relevant. Two of them are associated with the instability of hexagons to either the general solution (without rotation) 
or the oscillating hexagons (with rotation). We call the eigenvalues on these branches the 'Hopf eigenvalues.' The other 
two branches are associated with the translation modes of the hexagons; we call the eigenvalues on these branches 
'translation eigenvalues.' Note that they are identically zero only for (q,p) = (k, 0), i.e. for yj q 2 + p 2 — k and 8 = 0. 
For other perturbation wavevectors they are non-zero. In fact, in the presence of rotation the translation eigenvalues 
need not even be real but can form a complex pair. In the following we discuss the dependence of these four eigenvalues 
on the magnitude and, in particular, on the orientation of the wavevector of the perturbation. As discussed in the 
previous section, we restrict ourselves to the range —30° < 8 < 30° for the angle between the perturbation mode Pooo 
and mode A. 
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FIG. 4. Fourier space diagram with the shaded areas indicating perturbations Pooo to which the hexagons in the minimai 
model are unstable for R = 0.0294, k = 1.01, a = 0.1, 7 = 6. The sixfold symmetry is due to the rotation symmetry. In this 
case the instability first occurs at a finite angle 6. 

In Fourier space, the growth rates exhibit a sixfold symmetry of areas of instability, centered around the wavevectors 
of the hexagon modes (fc, 9 = 2irn/G), due to the rotation symmetry of the system (cf. p7j ). FigJ| shows the instability 
regions for Pooo in the vicinity of mode A for k = 1.01. In Fourier space (cf. Fig. [l]), a rotation of Pooo about A by 
60° is equivalent to a rotation of the system about the origin by 60° in the same direction. The clockwise rotation of 
Pooo about A by 60° makes its position equivalent to that of P010 before the rotation, the position of P010 after the 
rotation equivalent to that of Pj 10 before the rotation, and so on. Hence, the stability calculation exhibits a sixfold 
symmetry of unstable perturbation modes. Additionally, for arbitrary perturbation wavenumber (q,p), the real part 
of the translation eigenvalues is not necessarily zero when 8 = 0; this is only true when \J q 2 + p 2 = k. 

In the following we discuss the properties of the minimal and the extended model separately. In both cases we 
consider only those parameter values for which the hexagons are stable to perturbations within the hexagonal lattice. 
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This means that we consider values of the bifurcation parameter R between the saddle- node and the Hopf bifurcation 
points. Hence, the Hopf eigenvalues have negative real part for (q,p) = (k, 0) in all our calculations. 



A. Stability Properties of the Minimal Model 



In the minimal model (Q) we only found side-band instabilities involving the two branches of perturbations that are 
connected with the translation modes, i.e. the growth rates of the perturbations vanish as the perturbation wavevector 
(g, p) approaches that of a hexagon mode. For small R the instabilities are long-wave and the wavevector of the most 
unstable perturbation is close to that of a hexagon mode, i.e. the instabilities arise at a small angle 9 between the 
perturbation modes and the hexagon modes. The eigenvalues are real for small rotation rates, but become complex 
as the rotation rate is increased. 

Further away from threshold the instabilities become short-wave. This is shown in Fig.|| where the real part of 
the growth rate of the least stable modes is shown as a function of the angle 9 of the perturbation modes. In this 
illustration the magnitude of the perturbation wavevector is chosen to be equal to that of the hexagons, \/ q 2 + p 2 = k. 
For small angles the perturbations are damped. Only for larger angles the real part of the growth rate changes sign 
(dashed line). In this case the relevant eigenvalues are complex. It is worth mentioning that the translation mode 
and the Hopf mode (solid line) influence each other in this case for small angles. This will become more important in 
the extended model (cf. Fig.ff] below). 



CD 
CO 

rr 



O 



CD -o.oi - 




Angl 



5 

e of Perturbation 9 



FIG. 5. A short-wa ve instability driven by the translation modes in the minimal model for a = 0.1, 7 = 6, R = 0.04, and 
= 1.01 = y / q 2 +p 2 . 
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FIG. 6. Stability regions for the minimal model for various values of 7. Open circles denote sideband instability in all 
cases. Solid lines: neutral curve and saddle-node bifurcation. Dashed line: instability of hexagons to rolls. Dotted line: lower 
stability boundary of the rolls (cf. Fig. |^). a) 7 = 2. Solid symbols denote long- wave instability, b) 7 = 6. Crosses denote 
numerical instability of hexagons, diamonds denote stability limit of modulated oscillations, c) 7 = 14. Line denotes value of 
R of numerical investigation (see Fig.[l4| below) . 

The stability limits obtained from ( |l5| ) are shown in Fig.|| for a range of rotation rates 7. In these calculations 
the perturbation wavenumber is allowed to vary within the range 0.9 < \J q 2 + p 2 < 1.1 and the angle within the 
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range —30° < 9 < 30°. We note that variation of the non-rotational quadratic terms, in this case a, simply varies 
the strength of the subcriticality of the hexagon solution and shifts the instability of the hexagons to rolls up and 
down, thus setting the over- all scale for the stability range of the hexagons; no qualitative changes in the stability 
limits were found. The solid lines in Fig.^ indicate the neutral curve and the locus where the hexagons first arise in 
a saddle-node bifurcation. The dashed line indicates where the hexagons become unstable to rolls (via the unstable, 
general solution in which all three modes A, B, and C have different amplitudes (Fig. ||)). In the generic case in which 
the rotation affects the difference between the cubic coupling coefficients f% and this instability is oscillatory and 
leads to the oscillating hexagons. The dotted line indicates the lower stability limit of rolls with respect to hexagons. 
The side-band instability of the hexagons is shown with open circles. Its long-wave limit is given with small, solid 
symbols. As mentioned above, the long-wave result gives the correct stability limit for small R. For larger R, however, 
the short-wave instability sets in first. As the rotation rates is increased the stability limits become more asymmetrical 
as shown in Fig.pp. Also shown in this case (as crosses) are the stability limits as they are obtained in full numerical 
simulations of (0). For sufficiently large rotation rates the stability region of the hexagons detaches from the line that 
indicates the transition to the general solution. Thus, over a range of the control parameter R the hexagons are still 
stable with respect to rolls, but nevertheless at no wavenumber arc the hexagons stable with respect to side-band 
instabilities. This suggests the possibility of persistent dynamics, which is discussed below. 

B. Stability Properties of the Extended Model 

Within the extended model (|J) the coefficients ji and fs differ from each other for non-vanishing rotation rates 
and in the absence of side-band perturbations hexagons become unstable to oscillating hexagons with increasing R 
(dashed line in Figs. HJ§b). With side-band perturbations included it turns out that this Hopf bifurcation can occur 
first for finite angles between the perturbation mode and the hexagon mode, as shown in Fig.[?]a. Thus, depending 
on parameters the hexagons can become unstable with respect to modes rotated relative to the hexagons in two 
different ways: the destabilizing modes can be connected with the translation modes or the Hopf mode. There are 
two possibilities for the transition between the two regimes to occur. Either the two modes interchange their roles 
without affecting each other, or the transition involves a switching between the two branches. The latter is illustrated 
in Fig.0. While for R = 0.035 and k = 0.950 the destabilizing mode is clearly associated with the Hopf mode (Fig.^i), 
it is clearly connected with the translation branch for R = 0.020 and k = 0.957 (Fig.|^c). The transition occurs around 
R = 0.028 and k = 0.953 for which value the Hopf and the translation-mode branch merge and the destabilizing mode 
cannot uniquely be attributed to either branch (Fig.^b). For these parameter values the translation modes and the 
Hopf mode are in resonance, i.e. they have the same frequency. Fig|| shows where the stability limits are due to the 
Hopf branch and to the translation-mode branch, respectively. Comparing Figs.0a and ^c one may expect that the 
two types of instabilities would lead to different behavior. This is, however, not the case (see sec.|v| below). 
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FIG. 7. The branch switching transition in the extended model for a = 0, /3 = 0.2, 7 = 6, gi = 1 
k = 0.950 = y/q 2 +p 2 . An oscillatory instability driven by the Hopf modes, b) R = 0.028, k = 0.953 = ^fq 
and translation eigenvalues meet at two points; in this situation it is not clear wh ich dash -dotted eigenvalue branch belongs to 
the Hopf and which to the translation eigenvalues, c) R 
the translation modes (cf. Fig.B) . 



0.020, k = 0.957 = \J q 2 + p 2 . An oscillatory instability driven by 



In the absence of rotation, all eigenvalues are real. As rotation is increased, the translation eigenvalues become 
complex for a range of 9. In Fig.Ha a case is shown in which the translation eigenvalues are a complex conjugate 
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pair with negative real part for small \0\, but for larger \9\ the eigenvalues split and turn real with one eigenvalue 
becoming positive driving a steady instability. Here the only rotation term is g2- Thus, for 92 = there is no rotation 
in the system and all eigenvalues are real. One can get qualitatively similar instabilities with only 7 providing the 
rotation. However, rotation driven by 7 does not generally have a real eigenvalue becoming positive for such a large 
\6\ as shown in Fig.^Ja. Fig.|]b shows the resulting stability limits. As in the cases discussed above, close to threshold 
the limits are given by the long- wave perturbations (small filled symbols), while further above threshold the relevant 
instabilities become short-wave. In contrast to the cases discussed before, however, the short-wave instability now 
involves a single real eigenvalue. 
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FIG. 8. Stability regions for various values of g2 with a = 0, j3 = 0.2, 7 = 6. Solid lines: neutral stability and hexagon 
saddle-node bifurcation curves. Dashed line: Hopf bifurcation to oscillating hexagons. Solid symbols: instability associated 
with the Hopf mode, open symbols: instability associated with the translation mode. 
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FIG. 9. a) One of the real eigenvalues (thin dashed line) resulting from the splitting of the complex conjugate pair (thick 
dashed line) of translation eigenvalues is positive for R = 0.035, k = 0.938 = y/ q 2 + p 2 , a = 0, = 0.2, 7 = 0, (72 = 5. Note 
that individual eigenvalues can change from being real to being part of a complex conjugate pair of eigenvalues as 6 is varied. 
b) Stability regions for various values of gi with a — 0, (5 = 0.2, 7 = 0. Solid symbols: oscillatory instability; open symbols 
instability involving a real eigenvalue. 



IV. NUMERICAL SIMULATIONS 



To support the analytical calculation, and to study the nonlinear behavior ensuing from the instabilities, numerical 
simulations of (Q) and (||) were performed. A Runge-Kutta method with an integrating factor that computes the linear 
derivative terms exactly was used. Derivatives were computed in Fourier space, using a two-dimensional complex fast 
Fourier transform (FFT). 

The numerical simulations were done in a rectangular 32-7T by 647r/v3 box of aspect ratio 2/\/3 with periodic 
boundary conditions. The aspect ratio of 2/V3 was used to allow for regular hexagonal patterns. For the small 
amplitudes investigated here it was sufficient to use a 128 by 128 lattice keeping only up to the third harmonics of the 
modes on the critical circle. This allowed a larger number of modes within the critical annulus, reducing the effects 
of anisotropy. 

In all cases investigated, if the instability is associated with a real eigenvalue it does not saturate and eventually 
leads again to steady hexagons with a wavevector in the stable region. This is true for the usual long-wave case close 
to onset as well as for the case when the instability arises first at a finite modulation wavevector (cf. Fig.pk). 
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The oscillatory instabilities exhibit more interesting behavior. While they also can take the hexagon pattern back 
into the stable wavenumber band, over a range of parameters they lead in a supercritical bifurcation to modulated 
hexagons. They are quite similar to the oscillating hexagons identified earlier by Swift |J and Soward |ll| in that the 
hexagons are temporally modulated in time with modes corresponding to different orientations being phase-shifted 
with respect to each other. In the modulated hexagons the oscillations occur, however, predominantly in three side- 
band perturbation modes, which implies also a spatial modulation of the hexagon pattern. A typical time sequence 
over one period T = 225 is displayed in Fig.[n| Since the wavevector (q,p) of the sideband perturbation is close to that 
of the steady hexagons only one modulation wavelength fits into the system. Fig.[ll] shows that the spatio-temporal 
modulation is in the form of a standing rather than a traveling wave. The modulated hexagons are stable in a small 
region beyond the stability limit of the steady hexagons as delimited by the crosses and diamonds in Fig.^Jb and in 
Fig.0. 
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FIG. 10. A detail of stability diagram Fig.^|j. Circles: stability boundary of hexagons from weakly nonlinear analysis; crosses: 
numerical stability boundary of hexagons on the 128 by 128 lattice; solid diamonds: stability boundary of the modulated hexagon 
state. 
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FIG. 11. A real-space illustration of the stable oscillations of the perturbation mode for R — 0.0355, k — 1.02, a = 0.1, 
13 = 0, 7 = 6, <?2 = at times t = 0, t = 50, t = 110, t = 170. Indicated in the figure are the zero level contours of the solution. 

When the wavenumber k or the control parameter R is increased beyond the line marked by diamonds in Fig.^Jb 
the modulated hexagons become unstable and the regular pattern breaks down leading to a temporally and spatially 



chaotic state. Typical snapshots are displayed in Fig. 12i,b. They suggest a strong signature of a hexagonal symmetry 
and a predominant orientation of the disordered pattern which rotates in time. This is brought out more clearly in the 
power spectrum of Fig.|l2|a pictured in Fig.|l^c, which clearly shows six dominant peaks (see also Fig.|l^ below). Its 
temporal evolution is depicted in Figjl^a as a space-time diagram of the radially integrated power spectrum Fig.|l2|c 
with time increasing upward. Since ip is real, only the range —90° to 90° is shown. The rotation of the spectrum 
is seen to be surprisingly pronounced and steady over long periods of time, which are intermittently interrupted by 
periods during which the peaks in the spectrum are less pronounced. The space-time autocorrelation function of this 
state as obtained from a run of duration T max — 278, 880 is shown in Fig.|l3|b. The azimuthal decay of the correlation 
function is seen to be quite slow. We do not expect that the sixfold symmetry of the power spectrum persists as the 
system size is increased. In fact, the highest peaks in the power spectrum decrease by roughly a factor of two if the 
linear dimensions of the system are doubled. However, Fig.|l3|a,b indicates that the disordered state at least locally 
exhibits a structure with hexagonal symmetry, which rotates in time. 
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FIG. 12. a) Zero-level contours in real space for time t = 118580. b) Zero-level contours at a time when the power spectrum 
has rotated by about 30°. c) Power spectrum on a logarithmic scale of pattern at time t = 118580 for a — 0.1, 7 = 6, R = 0.036, 
k = 1.02. Note the residual sixfold symmetry. 




Angles Angles 
FIG. 13. Temporal evolution of radially integrated power spectrum and space-time correlation thereof in the disordered 

domain a = 0.1, 7 = 6, R = 0.036, k = 1.02 (cf. FigQ. 

The modulated hexagons and the disordered hexagons shown in Figjll] and Fig.[l2] are obtained in the minimal 
model (0) and arise from the short-wave instability connected with the translation mode. Very similar results are also 
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obtained in the more general model (|3j) where they can also arise from the short-wave Hopf mode. For the nonlinear 
evolution the long- wave behavior of the branch of destabilizing modes is apparently not relevant. 

An interesting aspect of the disordered state of FigJlJ is that it exists stably for the same value of the control 
parameter as the steady hexagons and the modulated hexagons. This coexistence of ordered and disordered states is 
reminiscent of the situation in Rayleigh-Benard convection with small Prandtl number where the state of spiral-defect 
chaos coexists with that of straight parallel rolls [2j. 

It is less surprising to find complex dynamics in parameter regimes in which no stable steady hexagons exist. Such 
a case is shown in FigJlJ. The parameters are as in Fig.^jc with R (marked by a line in Fig.^c) chosen in the gap 
between the stable steady hexagons and the transition line to unmodulated oscillating hexagons. In this case the 
disorder is sufficiently strong that the power spectrum is essentially isotropic for the system sizes investigated. 




From an experimental point of view an important question is whether the interesting dynamic states identified here 
are actually accessible even if the initial conditions cannot be prepared carefully. After all, they occur in parameter 
regimes in which the system could also evolve into a stripe pattern. We have thereforeperformed a few simulations 
starting from random initial conditions. For the parameters in the stability gap of Fig.||c the system evolves directly 
to a disordered state like that shown in Fig.|lJ. Fig.[l5| shows the space-time diagram for a simulation starting from 
small random initial conditions for parameter values for which steady and modulated hexagons coexist stably with 
the disordered state (cf. Figjl^). Clearly, the spectrum is initially essentially isotropic (bottom panel), reflecting the 
random initial conditions, and then evolves into a more ordered spectrum (top panel) corresponding to the rotating 
state discussed above. These runs suggest that the states found here should be experimentally accessible. Of course, so 
far no information is available to what extent the Swift-Hohenberg models discussed here describe a concrete physical 
system and whether the parameter regimes in which the dynamical states arise are accessible experimentally. 
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FIG. 15. Temporal evolution of radially integrated power spectrum starting from random initial conditions for a = 0.1, 
7 = 6, R = 0.03 (cf. Figjl|). 



V. CONCLUSION 



Motivated by the Kiippers-Lortz instability of roll patterns in the presence of rotation we have investigated the 
effect of rotation on the stability of hexagonal patterns. To get a first impression we have studied a minimal and 
a more general Swift-Hohenberg-type order-parameter model. In a linear stability analysis of the weakly nonlinear 
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hexagon pattern we focussed on instabilities that lead out of the space spanned by the modes corresponding to the 
hexagons themselves. Long- and short-wave instabilities were found, which can be steady or oscillatory. Immediately 
at onset the wavenumber band is limited by an instability arising from the translation mode. For larger rotation rates 
it may be oscillatory. In contrast to the case of roll patterns in the presence of rotation, which for sufficiently large 
rotation rates can become unstable at all wavenumbcrs immediately at onset due to the Kiippers-Lortz instability, no 
such parameter range is found for the hexagons. 

Further above threshold the instabilities typically become short-wave and oscillatory. They can lead to periodically 
oscillating spatially modulated hexagons or to a disordered state characterized by cellular and stripe-like elements. 
Locally the structure can exhibit a surprisingly regular rotation of the orientation of its residual hexagonal symme- 
try. It is worth comparing this behavior with that arising in the Kiippers-Lortz regime of systems with up-down 
symmetry. While in numerical simulations of an order-parameter model a very pronounced switching between the 
three orientations of the stripe pattern was found , in experiments also a behavior was observed that is more akin 
to the rotation of the spectrum found here ||. Particularly striking is the fact that the disordered state found here 
can arise for the same control parameter value as the periodically oscillating modulated hexagons and the steady 
hexagons. This apparent stable coexistence of the spatio-temporally chaotic state with an ordered state is reminiscent 
of the spiral-defect chaos found in Rayleigh-Benard convection, which coexists with the straight-roll state j34j . An 
interesting question is whether propagating fronts separating the ordered from the disordered state can be found as 
has been in those experiments. 

In investigations of the stability of hexagons with broken chiral symmetry within the framework of coupled Ginzburg- 
Landau equations the same type of long- and short-wave, steady and oscillatory instabilities have been identified 
recently |3f| . There also a regime was found in which the instability does not lead to an ordered state although such 
a state is stable for the same value of the control parameter. In the Ginzburg-Landau equations, however, that state 
cannot be described appropriately. They break the isotropy of the system and can only describe perturbations that 
make a small angle with respect to the hexagon pattern. Since the disordered state that was found in the Swift- 
Hohenberg equations investigated here exhibits an almost isotropic Fourier spectrum (cf. Figs.[l2|,[l4|), it leads to the 
excitation of ever higher Fourier modes in the numerical simulation of the Ginzburg-Landau equations. This renders 
them useless for the investigation of the chaotic state. 

In order to make concrete predictions for experimental systems the order-parameter equations have to be derived 
for specific physical systems. Situations of interest are Marangoni convection and Rayleigh-Benard convection with 
broken Boussinesq symmetry. For Marangoni convection with poorly conducting boundaries |^5|,[33| the effect of 
rotation has only been included recently in the derivation of long- wave equations Q . The case of Rayleigh-Benard 
convection has been treated to some extent. For poorly conducting boundaries a systematic long- wave expansion 
has been performed p4| . By imposing asymmetric boundary conditions (but neglecting other non-Boussinesq effects) 
quadratic terms in the order parameter equation were obtained, including the crucial term involving 7 that breaks 
the chiral symmetry. In that work the analysis focused on square rather than hexagonal patterns Jp4j . 

If the pattern arises first at a finite (i.e. not small) wavenumber an equation of the form (|^) or (gjcan be obtained 
in an approximation in which the kernel of the nonlinear integral term in Fourier space is approximated by local terms 
in real space |3(],|2(],|3l| . For Boussinesq Rayleigh-Benard convection with and without rotation it has been found that 
such an approximation can yield good agreement with weakly nonlinear amplitude equations while still preserving 
the isotropy of the system [Q. The non-Boussinesq case with rotation has not been treated yet. 

We gratefully acknowledge helpful discussions with B. Echebarria, A. Golovin, A. Mancho, W. Pesch, and M. Silbcr. 
The computations were done with a modification of a code by G. D. Granzow. This work was supported by D.O.E. 
Grant DE-FG02-G2ER14303 and NASA Grant NAG3-2113. 
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